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Abstract 

The characteristics of the electron-photon showers initiated by 2 to 10 GeV electrons 
aligned along the < 111 > - axis of tungsten crystals are compared with those for the 
amorphous tungsten . In this energy range, as known, the positron yield at the optimal 
target thicknesses is larger in a crystal case only by several percent. However, the 
amount of the energy deposition in a crystal turns out to be considerably (by 20 - 50%) 
lower than in an amorphous target providing the same positron yield, while the peak 
energy-deposition density is approximately of the same magnitude in the both cases. 
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1 Introduction 



The formation of electromagnetic showers in ahgned single crystals was actively studied dur- 
ing the last decade. The first experimental investigation of such showers has been performed 
in jl] at very high energy of incident electrons. Corresponding theoretical studies were started 
with ^ where an analytic solution of the problem was obtained, assuming that energies of 
all charged particles and photons involved are very high. This limitation was surmounted in 
[3J by suggesting a specific form of the radiation spectrum at axial alignment and performing 
corresponding simulations. Using this approach, the results of P for Ge crystal had been 
reproduced in^. The results of |2| are consistent with those of where another approach 
was used to obtain the radiation spectrum. In 0,0, and [U], the shower characteristics, 
such as spectral-angular distributions of photons and positrons as well as the amount of en- 
ergy deposition have been obtained depending on the kind and thickness of crystal targets. 
Investigations of the specific shower formation give good grounds for the idea proposed in 
[Zj, to substitute in a positron source an axially aligned crystal target for an amorphous 
one, as the enhancement of the radiation intensity is maximal just at the axial alignment. 
In further experiments (see [8-13]) using 1.2-10 GeV electrons aligned to the < 111 >- axis 
of tungsten crystals, measurements of some shower characteristics were already aimed to 
the development of a crystal-assisted positron source. Theoretical estimations performed in 
display a rather good agreement with results of recent experiments [10-13]. So, we can 
rely on our understanding of the physics of shower formation and on numerical results, at 
least for tungsten crystals in the energy range of incident electrons below 10 GeV. Note that 
just this energy range is proposed in future linear collider projects ( 2, 6.2, and 10 GeV 
correspondingly for CLIC [T3],NLC [12], and JLC [T7j ) and is considered here. 

Let us define the positron yield as the number of accepted positrons per one incident 
electron and the optimal target thickness as that providing the maximal yield. According 
to [3], [5], the maximal yield from a crystal target is always higher than that from an 
amorphous one and the excess goes up when the electron energy increases. However, the 
magnitude of such an enhancement is small, less than 14% even at 10 GeV. The more pro- 
nounced advantage of crystal targets appear in a considerable (by a factor of two at 10 GeV) 
decrease of the energy deposition. Indeed, the thermal effects caused by the energy deposited 
in a target are a critical issue for any powerful positron source based on the conventional 
scheme. We dwell mainly on this issue in the present paper. Below qualitative arguments are 
given explaining the lower energy deposition in crystals. The total deposited energy and the 
distribution of its density over the target volume are calculated for crystal and amorphous 
tungsten targets using the parameters of CLIC, NLC , and JLC . Thereby, a possible gain for 
these projects resulting from the use of crystal targets in the positron source is estimated. For 
accurate studies of thermal effects, some improvements have been performed in the computer 
code developed in 0, 0. The updated version of the code is used to study both crystal and 
amorphous cases. 

2 Energy deposition in crystal and amorphous targets 

In the energy range under consideration we are dealing with a "soft" (see j2]) shower when 
pair production is entirely due to the conventional Bethe-Heitler mechanism, while the crystal 
structure reveals in a considerable enhancement of the radiation intensity and a softness of 
the photon spectra. Remember that this enhancement decreases when the particle energy 
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does so as the shower develops. Starting with some depth Lq (see discussion in [O]), 
further development of the shower proceeds more or less in the same way for any (crystal or 
amorphous) type of the remaining part of the target. For the sake of simplicity, calculations 
are performed here for the all-crystal targets. However, they may serve as a good estimate for 
hybrid targets of the same total thickness and with a crystal-part length of the order of Lq. 
Let us remind that a hybrid target consists of a photon radiator made of a crystal followed 
by a pair converter made of an amorphous piece. From Fig^ a value of Lq > 0.2cm(0.57Xo) 
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Figure 1: Mean energy of the initial electron in units of the incident beam energy, Eq = 
10 GeV (a) and : fractions of the total energy carried correspondingly by all charged particles 
and photons (b). Solid lines are for amorphous and dashed for crystal tungsten target. The 
incident beam is directed along the < 111 >- axis of the crystal. 

can be chosen for 10-GeV electrons, since the fraction of the total energy carried by photons 
(~ 0.72) has been already saturated at this depth and the mean energy of the primary 
electron is sufficiently low to eliminate crystal effects in the last part. Such a saturation 
takes place in amorphous targets as well, but with a lower conversion level (~ 0.59) and at 
substantially larger depth (~ 5Lq at 10 GeV). Only a small part (less than 0.4% in the above 
example) of the beam power is deposited over Lq and the energy-deposition density is much 
less (about 8 times at 10 GeV) than its peak value. So, the crystal part of a hybrid target is 
not exposed to serious thermal effects which appear at larger depths in the later stage of the 
shower development. 

From calculations performed in [S], jB], the energy deposition in equivalent (providing the 
same positron yield) targets is always less in the crystal case. Let us present some qualitative 
arguments to explain this fact. The main process leading to the energy deposition is the 
ionization loss, the rate of which, q{z), reads approximately as q{z) ~ Cq ■ Nch{z), where 
Nch{z) is the number of charged particles at the depth z. Strictly speaking, the coefficient Cq 
may depend on z but its small variation as well as a small difference of Cg-values in crystal 
and amorphous cases are neglected in our estimation. So, the total energy, Q{L), deposited 
over the thickness L reads 

L L 

Q{L) = j dzq{z) ^Cgjdz N,n{z) , (1) 
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or, going over to the variable Nch{z) 



Q{L) ^ Cq 



dN. 



ch 



ch 



dz 



For sufficiently large L, the positron yield is roughly proportional to the total number of 
charged particles, Nch{L), i.e., for equivalent targets, the integrals in ^ are taken over the 
same region of variable Nch in both cases. To prove our statement, it remains only to verify 
that the logarithmic derivative which appears in the denominator of the integrand in Q is 
larger for crystals. This derivative, or logarithmic increment, characterizes the growth of the 
number of charged particles. As seen in Fig^ the conversion of the initial electron energy into 
photons is going faster for crystals, where, correspondingly, the pair-production process starts 
earlier and is more intensive, resulting in a larger increment. For the purposes of illustration, 
the energy-deposition rate per charged particle and the logarithmic increment are shown in 
Fig|21 We emphasize that the energy-deposition rate is practically the same for crystal and 
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Figure 2: Energy deposition rate per charged particle (a) and the logarithmic increment (b) 
in tungsten targets at Eo =10 GeV . Solid lines are for amorphous and dashed for crystal 
targets. 

amorphous samples, being almost independent of the electron energy. The non-constancy of 
this rate is mainly due to the contributions of other processes like photon absorbtion and 
annihilation of positrons which are taken into account by our simulations but were ignored 
in the estimation of q{z) used in ([Q). Evidently, the role of these processes increases with 
growing depth. A slow decrease of the rate at comparatively small depth is due to that of 
the mean particle energy. The point is that we use the so-called non-restricted energy loss 
description, where the rate diminishes when the particle energy does so. On the whole, the 
results of simulations presented in FigI21 con&m the above qualitative considerations. 

As shower develops, particles of sufficiently low energy may substantially change their 
direction of propagation due to single or multiple scattering. So that there even appears 
a "backward" flux of particles moving upstream; they do not increase the positron yield, 
however, they heat the target. This contribution is taken into account in our simulations, 
where we are able to trace separately effects from "forward" and "backward" particles and 
photons. To study the shower characteristics depending on depth (z), the target is divided 
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into slices by planes perpendicular to the incident beam direction, taken as 2;-axis. The 
spacing dz = 0.25 Xq is used (remember that Xq =3. 5mm for tungsten). Energy deposition is 
simulated within each slice and various distributions are recorded at the right-hand boundary 
( plane ) of the slice. In particular, the development of momentum distributions for positrons 
and photons is obtained as well as that of the beam spot size. Shown in Fig|21^a) is the energy, 
AEiiep{z), deposited in slices in units of Eq. This quantity is higher for the lower energy. 
Note that the deposited power, P^ep, reads as Pdep = Po • Edep/Eo where Pq is the incident 




Figure 3: Fraction of energy AEdep{z) / Eq deposited in ^-slices (a) and number of accepted 
positrons as a function of the deposited energy Edep/Eo (b) at Eq = 2 GeV ( curves 1 ), 
6.2 GeV ( curves 2 ),and 10 GeV ( curves 3 ). Solid lines are for amorphous and dashed for 
crystal targets. 

beam power. The contribution of the "backward" particles to E^ep increases with the target 
thickness L; it amounts typically to about 20% at L ~ AXq and thereby is not neghgible. The 
results of simulations presented in FiglHl (b) clearly confirm once more (cf. Fig.9 in the 
statement concerning the comparison between the values of Edep in the equivalent targets. 
The positron yield is calculated using the "theoretical" acceptance conditions from [12] ( 
see Figs. 11 and 12b in p!3^ ) in all three cases. The angular spread of the incident beam is 
neglected and the transverse size of the beam, a{ax = ay = a) is set to 1.6 mm at 2 and 6.2 
GeV, and to 2.5 mm at 10 GeV. At Eq=2 GeV, a simulation for (T=2.0 mm was performed 
as well. Corresponding results, not shown on graphs, are presented in the Table below. 

Let us remind now that, at equal depths and initial energies, charged particles are softer 
and have a larger angular spread for crystal targets (see, e.g.. Figs. 2-5 in J3] and corre- 
sponding discussion). All other things being equal, positron spectra are softer at a lower 
initial energy. As an illustration of these features, the mean energy and the transverse mo- 
mentum, < pt>i oi "forward" positrons are shown in Fig^a). We emphasize that, starting 
with L ~ Xq, < > is almost constant and practically independent of the initial electron 
energy and the type of the target (note the merging of six different curves on Fig^a)). Such 
result confirms and extends those concerning an amorphous target; it is essentially due to a 
counterbalance between the increase of the angle as a consequence of the multiple scattering 
and the decrease of the positron energy with increased thickness jS]. At the same time, the 
larger angular spread leads in the crystal case to a larger beam spot area, S", as seen in 
FiglH^b) where S/ Sq is plotted for "forward" charged particles. For evident reasons, the spot 
area of the "backward" charged flux ( not shown in FigE]) is somewhat larger (typically by 
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Figure 4: (a) : Mean positron energy (in MeV, upper curves) and transverse momentum (in 
MeV/c, lower curves); (b) : beam spot area development. Incoming energies are Eq = 2 GeV 
( curves 1 ), 6.2 GeV ( curves 2 ),and 10 GeV ( curves 3 ). Sq corresponds to the incident 
beam. The straight line y = 1 on (b) is drawn to guide the eye. Solid lines are for amorphous 
and dashed for crystal targets. 



20-25 %) than that of the "forward" one. 

Whereas the average deposited power can be handled somehow (e.g., by rotating the target 
and removing the excess heat through water cooling as used at SLC), the local and nearly 
instantaneous energy deposition is unavoidable, while being of critical concern for target 
damage. Basing on the analysis of the SLC damaged target ( see jS]), it is now adopted 
that the peak energy-deposition density should not exceed 35 J/g to ensure a sufficiently 
long term of safe operation. In our simulations, the total scanned volume is a cylinder 
coaxial with the incident beam direction. The radius, R, of the cylinder was about R ~ 
3.7(7, in which case less than one thousandth of initial electrons does not hit the cylinder at 
entry. This cylinder is sliced into disks of the thickness dz = 0.25 Xq, i.e., 2;-spacing is the 
same as in above calculations, allowing a mutual checking. In turn, each disk is divided by 
circles of uniformly increasing radii with the step dr = 0.02R into 50 parts - one internal 
disk (altogether, such disks form the internal cylinder ) and 49 rings. More precisely, we 
have (ir(cr=1.6mm)=0.12mm, dr {a =2. 0mm) = 0.15 mm, and (ir((T=2.5mm)=0. 185mm. The 
energy deposition and the number of charged particles is simulated in each meshed volume 
providing corresponding transverse distributions for each z-slice. It is noteworthy that z- 
dependencies of shower characteristics derived from 3-dimensional distributions coincide with 
those obtained in direct calculations, thereby verifying the consistency of the two essentially 
different forms of the output which we have used. From Figj^ (b), the energy- deposit ion 
density (EDD) drops at large cell radius, being maximal in the internal disk of each slice. 
Shown in Fig|Sl (a) is the EDD per one incident electron in the internal cylinder, versus the 
depth. Like the positron yield and the total number of charged particles , the peak values of 
EDD (PEDD) in crystals are somewhat higher and are reached at smaller depths, compared 
with amorphous targets. At equal beam size, PEDD is higher for the higher energy (cf. 
pairs of curves in Figl^l^a) calculated at 2 and 6.2 GeV for the same a =1.6 mm). However, 
by increasing the beam size, lower values of PEDD may be obtained even at higher energy. 
This is seen in Figl^ (a) if we compare the curves calculated at 10 GeV ( a =2.5 mm) with 
those at 6.2 GeV. Examples of transverse distributions of EDD are plotted in FiglSt^b) for 
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Figure 5: (a): Energy- deposit ion density in the internal cylinder, versus the depth, at 
indicated initial energies; (b): Transverse distribution of the energy-deposition density at 
z = 4.0X0. Solid lines are for amorphous and dashed for crystal targets. 

the same depth z = AXq. The distribution of charged particle density (not shown) is similar 
to that of EDD with a ratio EDD/density being almost constant within the first ten cells 
and roughly equal to the energy-deposition rate. In Figl^l^b) this similarity reveals in larger 
width of distributions for crystal targets (cf. FigEfb) ). The EDD is presented in Figl^in 
units of GeV ■ cm~^ which for tungsten corresponds to 1 GeV ■ cm~^ = 8.30 • 10~^^J/g. 

Due to an extremely short duration of the pulse, contributions to EDD from all the 
incident electrons are added and the resulting EDD-value is simply the product of the ob- 
tained EDD-value per e~ and the number of electrons per pulse, provided that the target 
reverts to initial thermal conditions during the repetition period. Some results concerning 

Table 1: Energy deposition in crystal and amorphous targets. The fraction of an 
incident beam power deposited in a target, Rth = Pdep/Po'i the peak energy deposition density 
,PEDD, in units of J/g; the energy,ii^o; and transverse size of an incident beam,cr; target 
thickness, L, measured in conventional radiation lengths, Xq; gain in the total deposited 
power, G = (1 — E^^^/E^^) ■ 100, from comparison of Wcr and Wam targets 



Beam 


Eo =2 GeV,a =2.0 mm 


Eo =6.2 GeV,(T =1.6 mm 


Eo =10 GeV,a =2.5 mm 


Target 


W75Re25 














Wam 


W,r 


L 


4.0 


4.0 


3.0 


4.0 


4.0 


2.5 


6.0 


6.0 


3.0 


Rth 


0.248 


0.238 


0.193 


0.142 


0.147 


0.106 


0.310 


0.291 


0.137 


PEDD 


35.0 


32.3 


33.8 


35.0 


56.1 


60.3 


35.0 


28.2 


30.6 


G 


19% 


28% 


53% 



the energy deposition in crystal and amorphous targets are presented in Table H For the 
sake of comparison, the thicknesses of amorphous tungsten targets (Wam) are the same as 
proposed for Wj^Re2^ targets in the projects [T3], ||16j, and |17j. For crystal targets iWcr), 
the corresponding thicknesses are determined from the equivalence rule (the same positron 
yield as for Wam target). Our values for PEDD given in Table [T] are obtained using the 
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same numbers of electrons per pulse, Ne{2 GeV) = 2.08 ■ 10^^ N^{6.2 GeV) = 0.96 • 10^^ and 
Ne{lQGeV) = 0.64 • 10^^ as in [HI, [IHl, and [TTj respectively. 

Since the parameters of the two amorphous targets such as radiation lengths and densities 
almost coincide, the same is expected for shower characteristics at equal depths. However, 
one should bear in mind that two different computer codes were used in simulations, EGS4 
for Wj^Re2^ targets, and our code for Wam targets. From Table HJ the concordance of the 
two approaches is almost perfect for the fraction of the deposited power, Rth- Concerning 
the peak energy-deposition density, PEDD, the values obtained for the Wam target at 2 and 
10 GeV are somewhat smaller than those for the W75-Re25 target and the distinctions are not 
too big. On the contrary, our value for the Wam target is substantially larger at 6.2 GeV. 
Let us argue the point in detail. At equal beam size, the peak EDD/electron is expected to 
be roughly proportional to the initial electron energy, Eq. This assertion follows from the 
qualitative consideration performed above and is verified by calculations. For example, if we 
compare the curves in Figl^fa) calculated at 2 and 6.2 GeV for the same a =1.6 mm, we 
obtain K =2.7 as the peak EDD ratio instead of K =3.1 from the rough estimate where 
K{Eq2, Eqi) = Eq2/Eqi. At L = 4Xo, which corresponds to NLC conditions, this ratio 
further diminishes up to 2.55 since the peak EDD value is achieved for Eq =6.2 GeV at a 
larger depth. So, at equal beam size, the relationship between PEDD values at different 
energies reads roughly as 

PEDD{2) ~ PEDD{1) ■ K{Eo2, ^oi) • iVe(^o2)/iVe(^oi) ■ (3) 

The PEDD in Wi^Re25 target at Eq =2 GeV was calculated in JHj not only for a =2.0 mm, 
but also for a =1.6 mm, where the value of 53.1 J/g has been obtained ( at these conditions 
we have 47.7 J/g for the Wam, target). Using the estimate Q with PEDD{l)=b'i.l J/g and 
K =2.55, we obtain for NLC conditions PEDD{2)=62.5 J/g which is consistent with our 
result for the Wam target. Conversely, starting with results of [T7j for the Wj^Re2^ target, 
we obtain the estimate {PEDD{l)='ib J/g, Em= 10 GeV, K =0.62) PEDD{2) ~ 33 J/g at 
a =2.5 mm, i?02=6.2 GeV. So, the PEDD value of 35 J/g at Eq=Q.2 GeV is more consistent 
with the beam size of a =2.5 mm than with a =1.6 mm indicated in pT^ . 

Comparing the magnitude of thermal effects in Wam and Wcr targets providing the same 
positron yield, we conclude that , using crystal targets, the total deposited energy can be 
considerably diminished while the peak value of the energy-deposition density is kept ap- 
proximately on the same level as in the amorphous case. 
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